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ABSTRACT 


In this thesis, mathematical models are established for 
the development of oxygen toxicity in divers. The study 
endeavors to derive the shape of the oxygen tolerance curve 
in terms of depth-time limitations by statistical analysis 
of existing data. By assuming a known distribution for 
the time-to-serious~-symptom, mathematically predictive 
models are developed which allow a greater degree of pre- 
dictability in mission profiles and allow the associated 


risk to divers to be evaluated. 
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IT. 2RBRODUCTION 


Priestly, upon discovering oxygen in 1775 philosophized 
whether, "it might exhaust the animal powers too soon, and 
as a candle we may burn too fast if used in our normal state 
of health" [Ref. 1]. It is known today that oxygen consump- 
tion in animals and man differs from Priestly's burning candle 
in that increases in tissue oxygen pressure do not accelerate 
Oxidative metabolism. However, Priestly's suspicion, that 
oxygen might have a toxic or deleterious effect, has proven 
wo be correct. 

The toxic effects of oxygen on the central nervous system 
(CNS) have been recognized since the later half of the 19th 
century. Paul Bert is given credit for the earliest investi- 
gations in this area [Ref. 2]. In his work "La Pression 
Barometrique" he described in detail the incidence of convul- 
Sions in various animal species exposed to high pressures of 
oxygen. In 1878 he subjected larks to 15 to 20 atmospheres 
of air and noted that they first convulsed and then died. 
Further studies by Bert showed that it was the oxygen rather 
than the nitrogen in air that caused this phenomenon. Lorrain 
Smith [Ref. 3] is another early pioneer (1899) in the field 
of oxygen toxicity. He carried out the first extensive inves- 
tigations of the toxic effects of oxygen on the lungs (pul- 


monary pathology). 





Since these early investigations, oxygen has been shown 
to have toxic effects on many tissues in the body when present 
in abnormally high partial pressures. For a detailed account 
of the studies up to the early 1940's, the reader is referred 
to Stadie, Riggs and Haugaard 1944 [Ref. 4] and Bean 1945 
(Mee. 5). 

During the last few decades, research in the area of oxy- 
gen toxicity has reached a high level of activity in response 
to increased applications in medical therapy and in diving 
operations. The use of oxygen to decrease the duration and 
increase the effectiveness of decompression procedures in mili- 
tary and commercial diving operations has resulted in the 
exposure of a large number of healthy men to high pressures 
of oxygen, seemingly without any harmful consequences. More- 
over, with research and experimentation a better understanding 
has evolved of the universal dependence of vital biological 
processes on oxygen, and on the opposing cellular antioxidant 
defense mechanism that balance the toxic effects of oxygen 
[Ref. 6]. It is now known that the same oxygen pressure re- 
quired to sustain life would cause lethal oxygen poisoning in 
the absence of these defense mechanisms [Ref. 7]. 

The severity of oxygen poisoning appears to increase pro- 
gressively with the increase in the partial pressure of oxygen 
and with greater exposure time. At sufficient pressure and 
exposure duration, oxygen will cause initial functional im- 


pairment and ultimate destruction of any living cell [Ref. 8]. 








The precise mechanism through which oxygen exerts its detri- 


mental effects on the tissues and organs of the body is the 


subject of continuing research and controversy. Among the 


numerous biochemical and biophysical explanations which inves- 


tigators have identified as possible causes are: 


er ) 


(=) 


Inactivation of enzymes, especially those containing 
sulfhydryl groups. 

The formation of powerful oxidizing radicals. 
Reduction of the amount of gamma amino butyric acid 
(GABA), a transmitter at CNS inhibitory synapses. 


Disruption of cellular membrane functions. 


The United States Navy's interest in the effects of hyper- 


baric oxygen stems largely from the following applications 


of oxygen in the conduct of diving operations: 


el) 


(2) 


As a breathing medium employed alone in the conduct 

of closed-circuit clandestine SCUBA operations. 

When used in this manner, oxygen has several distinct 
advantages over compressed air--it leaves no trail 

of bubbles by which a diver could be detected, and the 
absence of inert gas allows all the gas carried by the 
diver in his cylinder to be utilized, thus prolonging 
the operating time available. 

As a component of the atmosphere in mixed gas breath- 
ing systems that enable the Navy's divers to operate 
at deeper depths and for longer periods of time than 


1s possible with compressed air diving. 





(3) As an adjunct to decompression when used by divers 
who have incurred a decompression obligation breathing 
air, helium-oxygen or other mixtures. The divers are 
able to shorten their decompression time by breathing 
100% oxygen at their decompression stops and thus 
increase inert gas removal from body tissues. 

(4) As a component of the breathing mixture in saturation 
diving where oxygen is employed at higher than normal 
partial pressures to decrease the adverse effect of 
inert gas and to shorten decompression time. 

(5) Hyperbaric oxygen used in the treatment of decompression 
sickness and air embolism, either as prescribed in 
the standard Treatment Tables 5 and 6 or as part of 
Pietacimducilit zed Saturation treatment used for cases 
Vmehetall to cespend to conventional therapy (U.S. 
Navy Diving Manual) [Ref. 9]. 

The guidelines for oxygen diving in the U.S. Navy are 
published in the U.S. Navy Diving Manual. The maximum depth 
for closed circuit pure oxygen SCUBA diving during normal 
operations is 25 feet. The current depths and times for 
100% oxygen are shown in Table l. 

If one studies the existing data or Table 1, one observes 
that as the depth of the oxygen exposure increases, the dura- 
tion of allowed exposure decreases. This relationship 
reflects the fact that CNS oxygen toxicity is directly related 


to the partial pressure at which oxygen is breathed. A list 





DABBLE. 1 


OAYGEN DEPTH-TIME LIMITS 


1. Normal Operations - Depth (Ft) Time (Min) 
10 240 
15 Lae 
20 ILC 
2 fis. 
2. Exceptional Operations - Depth (Ft) Time (Min) 
30 45 
5S 25 
40 10 


Source: U.S. Navy Diving Manual 


of the factors which are known to affect the onset of oxygen 
toxicity 1S presented below. This list and much of the back- 
ground information on experiments was provided by Dr. Frank 
K. Butler of the Navy Experimental Diving Unit (NEDU) in 
Panama City, Florida: 
(1) Partial pressure of oxygen--the greater the depth 
Or pressure, the greater the likelihood of developing 
oxygen toxicity. 
(2) Duration of exposure--increases toxicity with increasing 
time of exposure. 
(3) Individual variations in susceptibility--this variation 


applies not only to the difference in susceptibility 
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between individuals but also for differences in the 
Same individual at different times. 

(4) Immersion in water--increases the onset of toxicity. 
Divers being treated for air embolism or decompression 
Sickness routinely breathe oxygen in a recompression 
chamber at 60 feet. Breathing oxygen at a depth of 
60 feet in the water would result in the fairly rapid 
onset of CNS oxygen toxicity. 

(5) Hypercarbia--divers breathing an increased partial 
pressure of carbon dioxide appear to be much more 
Susceptible to CNS oxygen toxicity. 

(6) Exercise--even in the absence of an exercise-induced 
hypercarbia, the exercising diver is more likely to 
incur toxic symptoms than the diver at rest. 

The most dramatic sign of oxygen toxicity is the grand 
mal convulsion. Convulsions may occur after other signs or 
symptoms have been experienced or they may be the first indi- 
cation of oxygen toxicity. If the hyperbaric exposure to oxy- 
gen is stopped, the convulsions subside after a variable period 
of time and are not thought to be associated with any permanent 
neurological damage. The diver must however, be protected 
against physical trauma secondary to the convulsions and 
possible air embolism as a result of being brought to the 
surface while unable to carry on normal breathing. And of 
course, any military mission in which this would occur would 


be jeopardized. A variety of other CNS symptoms may be 
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encountered. The mnemonic acronym suggested by the Navy 
Diving Manual is "Ventid": 

VISION. Visual abnormalities, especially "tunnel vision", 
which is a decrease in the peripheral field of 
V isaen:. 

mee. Audatory distwebances;, espec#ally tinnitus whch 
refers to a sound perceived by the ears but not 
Originating from an external source. The sound 
may be described as a ringing, buzzing, or 
machinery type of sound. 

NAUSEA. 

TWITCHING. The facial muscles are especially prone to 
Vapeching bub Ltemay occur in other groups of 
muscles. 

IRRITABILITY. A change in mental status which may be 
experienced by the diver as anxiety, confusion, 
Saeinabwiwry toeteecus his thoughts. 

DIZZINESS. 

meelo42eand 1943, a British investigator, H.K. Donald, 
conducted over 2000 man-dives for the purpose of establishing 
guidelines regarding oxygen toxicity for British frogmen in 
World War II. His research remains the largest dive series 
conducted to date investigating CNS oxygen toxicity. Some of 
his results and those of other earlier experimenters will be 


recounted. 
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Donald immersed 100 volunteers suited in Davis Submarine 
escape apparatus in a wet chamber at 50 feet of sea water 
(FSW) for a maximum time of 30 minutes. 30 FSW is about two 
atmospheres absolute (ATA). The water temperature was about 
65 degrees F. Of these 100 divers, 26 convulsed, 24 had other 
Symptoms of oxygen toxicity and 50 had no signs or symptoms. 
Donald also attempted to find the depth at or about which con- 
vulsions would not be observed. He conducted many such 
Similar experiments at varing depths [Ref. 10]. 

An NEDU study was published in 1947 by Yarborough, Behnke 
et al [{Ref. ll]. They tested various depths ranging from 30 
to 100 FSW. Ina resting dive underwater at 60 feet, 32 out 
eee 0/7 Cxposures were terminated prior to 60 minutes, which 
was the maximum time for the dive. The average time to 
termination was 32 minutes (ranging from 8 to 58 minutes). 

Two convulsions were noted: one at 13 minutes and the other 
at 24 minutes. Another interesting observation from this 
study was that subjects immersed in aay to their necks 
tolerated oxygen for periods similar to those of subjects in 
a dry chamber. 

Another study was conducted at NEDU in 1953 by Lamphier 
[Ref. 12]. He used a continuous flow oxygen facemask apparatus 
in 80 degree F. water. Work rate was set at "greater than a 
man would voluntarily sustain under diving conditions, although 
less than the maximum possible." No convulsions were observed 


at depths less than 35 feet. One convulsion occurred at 35 
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feet after 42 minutes (out of a total of 5 dives to a maximum 
time of 43 minutes). One convulsion was also noted at 40 
feet after 31 minutes (total of 13 dives to a maximum time of 
30 minutes). This convulsion at 3l minutes is as reported 

by Lanphier indicating that at least one dive went beyond 

30 minutes maximum. No convulsions were noted in 5 dives at 
45 feet for a maximum time of 19 minutes. 

In 1980 an NEDU study was published by Piantadosi, Clinton 
and Thalmann [Ref. 13]. Twenty-four oxygen exposures lasting 
from 80 to 271 minutes were performed by six immersed exer- 
Ccising divers at 25 FSW (1.76 ATA) in both warm and cold 
water. In this experiment two types of work were performed, 
moderate work (50 watts) for long periods of time, and graded 
(25-150 watts) lasting 85 minutes. CNS oxygen toxicity was 
not observed during these experiments although two divers 
had clinical evidence of early pulmonary oxygen toxicity. 

The latest series of oxygen studies conducted at NEDU were 
directed by Butler and Thalman (1982-1983). These experiments 
were conducted in the NEDU Ocean Simulation Facility under 
rigorous experimental conditions. Divers were exercised on a 
horizontal bicycle ergometer placed on a platform four to five 
feet below the water level. A series of variable depth dive 
profiles were conducted to probe the tolerance curve for 
reasonable profiles which allow for moderate depth excursions. 
The results of these experiments were not published as of the 
writing of this thesis but will be available through the Navy 


Experimental Diving Unit, Panama City, Florida. 
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In order to make the best use of oxygen, it is clearly 
necessary to have accurate information about safe durations 
of exposure to various partial pressures. Rules such as the 
Old 30-foot for 30-minute limitation on the use of 1003 oxy- 
gen in closed circuit diving are not very useful even if 
perfectly valid. Such rules give no information about the 
possibility of spending a shorter time at a deeper depth or 
a longer time at a shallower depth or the probability of an 
occurrence of a serious symptom which could abort a sensitive 
military diving operation. The basic need is limit curves of 
partial pressure versus permissible exposure times, showing 
the limits which apply to specific diving conditions for the 
degrees of safety desired expressed as probabilities. For 
this, the phenomena of oxygen toxicity must be mathematically 
modeled. 

Although the duration of exposure required to produce CNS 
oxygen toxicity is basically a function of the partial pres- 
sure of oxygen in the atmosphere--the higher the partial pres- 
sure the less time 1s required to produce symptoms--other 
factors are also very important. While the basic relationship 
of partial pressure and time hold for any specific situation, 
it could be overshadowed by the importance of these other 
variables which can greatly modify the tolerance curve. Of 
these, physical exertion and excess carbon dioxide appear to 
be the most important. Temperature, anxiety and a number of 


physiological factors might also be significant, but less 
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Strikingly so. The situation is further complicated by the 
fact that individual variations are large and almost entirely 
unpredictable [Ref. 10]. Not only do individuals differ from 
each other, but their tolerance changes from dive to dive 

for reasons which are unknown. This factor in concert with 
the lack of control in the experiments which have been con- 
Gucted thus far result in confusion and difficulties in 
handling the existing recorded data. 

Studies conducted in the past have established the signi- 
ficance of factors such as exertion [Ref. 10] but have not 
provided sufficient quantitative data to permit establishing 
a tolerance limit under the specified condition with any 
@egqree of confidence. In setting such a limit, all of the 
relevant conditions likely to be encountered deserve considera- 
tion. Individual variations in tolerance must be taken into 
account. The limits must be very safe because oxygen convul- 
Sions in SCUBA diving could easily end in death, or abort a 
sensitive military diving operation. However, making it 100% 
safe for all individual and all conditions is too restrictive 
and is out of the question. 

An ideal study of oxygen tolerance would include investi- 
gation of all significant variables, their degrees, and their 
combinations. A very large number of subjects would be em- 
ployed to illuminate the question of individual variability, 
and all exposures would be carried to convulsion, which is 


the only unequivocal end-point of oxygen toxicity. 
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An ideal study would also include determination of oxygen 
tolerance using a variety of gas mixtures at appropriate depths, 
Since the toxicity of oxygen may be modified by the presence 
of other gases. Studies have indicated that increased nitrogen 
pressure can suppress the convulsive manifestations of oxygen 
toxicity in animals. However, later NEDU studies have indi- 
cated that the presence of nitrogen may decrease oxygen 


tolerance in man [Ref. 12]. 


A. NATURE, SOURCES AND AVAILABILITY OF DATA 

This investigative effort was undertaken to explore the 
.feasibility of developing mathematical models which character- 
1ze the oxygen tolerance of human beings. Therefore, a survey 
of the scientific literature was performed to locate experi- 
mental data on time limits for human exposures to various 
partial pressures of oxygen. Unfortunately, data on human 
experiments published in the literature is extremely variable 
and 1s subject to influences, both physiological and environ- 
mental. Each scientist used a different approach to his experi- 
ments. Different breathing apparatus were used, different 
temperatures of the water, different work loads and depths, 
and other variants were explored by each scientist. This fact 
makes it extremely difficult and risky to compare the findings 
of the various experimenters. 

Raw data from the early experiments which were conducted 


by the Royal Navy Experimental Diving Unit during World War 


L7 





II is very difficult to acquire and some important experimental 
runs are only available as summaries of the results. A sum- 
mary of the findings such as Reference 10 is readily available 
but the raw data is unavailable, being lost or buried ina 
peund Of information in London, England. 

Most of the raw data on U.S. Navy experiments is available 
through the NEDU library in Panama City, Florida. This data 
is in the form of the actual log books and is generally avail- 
able in micro-fische. The result is the raw data is extremely 
tedious to work with and is often sketchy or unreadable. 
Fortunately, good summaries of the raw data are available in 
the form of NEDU reports such as References ll, 12, and 13. 

During the last few decades, with the increased interest 
in oxygen toxicity, a great deal of information has become 
available. Most of this information is of a theoretical or 
clinical nature and is not applicable to statistical analysis. 
It is very time consuming but informative to sort out the 
theoretical data from the empirical data and it is obvious 
that improved human data on oxygen exposure limits is sorely 
needed if statistical analysis is to develop a model with any 


high degree of confidence. 
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IIT. MATHEMATICAL MODELS 


In making a mathematical model for a real-world phenomenon 
it 1s necessary to make certain simplifying assumptions so 
as to render the mathematics tractable. Many different types 
of mathematical functions can be used to model a response 
that is a function of one or more independent variables. A 
mathematical model is merely a mathematical representation of 
a phenomenon or process. Ideally, we want a model to incor- 
porate everything that is important and significant about the 
phenomenon and ignore everything else. This buys the sim- 
plicity needed to provide a convenient method for performing 
certain computations. 

On the other hand, we cannot make too many simplifying 
assumptions, for then our conclusion or predictions obtained 
from the mathematical model would not be applicable to the 
real-world phenomena. A mathematical model is an abstraction 
which associates parameters and processes of the real-world 
with expressions and operations in a mathematical structure. 
Thus, we must make enough simplifying assumptions to enable 
us to handle the mathematics but not so many that the mathemati- 
cal model is no longer tractable. 

Keeping in mind that a density function is a theoretical 
model for a population of real data that occurs in the real- 


world, how do we know which model to use, and to what extent 


Jes, 





does it matter if we have an erroneous model? Any fitted 
model we could develop is merely an approximation, either 
because some variables are not recorded or incorrectly 
measured, or because the functional form used is not exactly 
correct. The mathematical model is merely a useful fiction 
suggested so that straightforward techniques can be used for 
the description of a relationship or the prediction of future 
values. It is unlikely that one would find a density function 
that provides a perfect fit, which would represent a perfect 
fit of the density function to nature. The purpose of a 
probabilistic model is to provide the mechanism for making 
inferences about the real-world; therefore a good model is one 
that yields good inferences about the real-world phenomena of 


interest. 


Peete DEL FOR THE DISTRIBUTION OF TIME TO SYMPTOM 

A reasonable selection of a model is sometimes implied by 
theoretical considerations. A second way to select a model is 
to form frequency histograms or empirical distributions of the 
raw data and choose a density function that would visually 
appear to give a similar frequency curve. One simplifying 
assumption that mathematicians are particularly fond of and 
that is often made is to assume that the random variables are 
exponentially distributed. The reason for this is that the 
exponential distribution is both relatively easy to work with 


and is often a very good approximation to the actual distribution. 
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Figure 1 shows the empirical density function for Donald's 
experiments [Ref. 10], at 103, 93, 83 and 73 FSW as computed 


by an APL function which smoothes the raw data into an empiri- 


cal distribution [(Ref. 14]. These experiments were chosen 
because they were carried to deciSive end-points. These 
empirical density functions look approximately gamma. If one 


truncates these graphs on the left to account for the fact 
that some time is required to press the divers to depth, the 
time-to-symptom is plausibly exponential. For those non- 
divers, the time of the dive begins and is calculated from the 
time the diver leaves the surface and it takes some variable 
amount of time depending on the individual diver to descend to 
the desired depth. Even with this limited number of pies, points 
there is a conspicuous build-up of hits followed by a long 
drawn-out tail. These graphs suggest that the exponential 
distribution could be used as a reasonable approximation to 
the real-world phenomenon of oxygen toxicity. Future analysis 
could explore the possibility of fitting gamma distributions 
to the times-to-symptom at depth and attempt to integrate 
these distributions in terms of a simple formula. 

The approach taken in this thesis will be to analyze the 
available data with existing mathematical models. For this 
effort certain mathematical tools will be necessary. For- 
tunately, many of the tools used in systems test and evalua- 
tion and systems reliability are applicable. 

In the study of the reliability of systems it is customary 


to think of the time-to-failure, up-time or lifetime of an item 


on: 
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as a random variable, T. Fundamental to this concept is the 
idea that the time-to-failure varies randomly from individual 
item to item and failure is a well defined notion. As an 
example, failure in a human lifetime could be defined as 
death which, with respect to diving, 1S an unambiguous con- 
cept. In oxygen toxicity studies failure could be defined as 
the time-to-convulsion or the time-to-symptom. In this thesis 
we will work with the time-to-serious-symptom. This is be- 
cause serious symptoms have been more accurately recorded and 
such toxic episodes usually resulted in the experimenter 
terminating the experiment. 

Let T equal the time-to-serious-symptom, then given that 
T 1s a random variable, reliability has the following 


definition: 
RiGee (tT > tt) 


Mien 1S, R(t) is the probability that T, which is the time-to- 
serious-symptom, is greater than t, where t is the duration of 
a particular misSion or mission time. If such a phenomenon 
is stable then such a measure has meaning in terms of a dis- 


meeoieclon function of T, F(t) : 


R(t) = P(T>t) = Ll- F(t) = Fim ve) 


Bete) Po—ereomemeconsiderca aswa £LUnction of the mission time, 
t, and is called the survivor function. Accordingly, mission 


reliability is the probability that a diver will not have a 
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serious symptom or convulsion during a specified period of 
time under specified environmental conditions. 

According to the above definition of mission reliability, 
the probability that a diver will not incur a serious symptom 
Or convulsion which would cause the abortion of a diving opera- 
tion, depends on the length of time a diver is at depth on pure 
oxygen. Thus of fundamental importance is the distribution 
of the times to failure (serious symptom or convulsion). A 
useful way to characterize this distribution is by means of 
its associated hazard rate of instantaneous failure rate. 

Suppose that a distribution function F(t) exists and 1s 


stable over time. This implies that: 


Now suppose that diver X has survived (no serious symptom or 
convulsion) for some time T, and we wish to know the proba- 
bility that he will not survive for an additional time dt, 
i.e., he will have a serious symptom or convulsion in the 
Subsequent small time interval dt. Let A(t) equal this hazard 


rate. Then: 
Mee) ser (De(t,t+at) |T >t) 


This says that A(t) is the probability that T is an element 
Sec ne Mmiterval trom & to — plus dt, given that T is greater 


Eman ct. By @detinition: 
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A (6) 


Pie e (eet et) ,T > t) 





a A(t) 


rl Sc 
. f(t) 
F(t) 
: -d/dt F(t) 
F(t) 
_ ~ (F(t) ]' 
Fn (t) 
; -~d(ln eae) ) 
— ae 
tC —— 
5 4! A(t)dt = In F(t) 
e 
=} eGe) che 
F(t) = e 
Wow R(t) = F(t) 
where /A(t) 


the distribution of failure times is: 


f(t) 


= 
= |) Mee) Ge is the integrated hazard rate, and 
0 


ie 


= Ce jicie 


= A(t)e : 


Since f(t) = )\ (eePR(t) 
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From the formula — 
the solution to this 


tt f£ollews that: 


du 


asin): 


differential equation is: 


Solving for ln F(t) by 
integrating both’ sides: 


Take the exponential of 


both sides 
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We know that the survival of diver X is a function of the 
partial pressure of oxygen (PPO2), duration of exposure, 
metabolic rate, partial pressure of carbon dioxide (PPCO2), 
anxiety, and other variables. Thus the integrated hazard 


could be written as: 


A(t) = A(VyiVor-eerVy) 


Mission reliability can then be written as: 


-A(v 


Ve) 
Rae = e a 


LiYorecee 
The probability that diver X gets a hit (here hit means the 
occurrence of a serious symptom or convulsion) is 1-R(t). 

The question is what is A(VyrVor-eeerV,)? Since an ideal study 
of oxygen toxicity has not been conducted, one must deal with 
those variables for which data is available namely, duration 
of exposure and depth, or depth as expressed as PPO2. To some 
extent metabolic rate has been explored in working versus non- 
working dives. These variables are of limited value due to 
the lack of commonality in control of different experiments 
and by the fact that PPO2, PPCO2 and duration of exposure at 


depth were not precisely measured in past experiments. 
B. EXPECTATION AND THE PROBLEM OF CENSORING 


The mean-time-to-failure (MTTF) iS a quantitative assess- 


ment of the failure times. 
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MTTF 1s an estimator of the expectation of the mean-time-to- 
failure. Inoxygen toxicity studies there are problems in 
calculating the expected value of T, E(T), which result from 
the termination of experiments prior to all the divers reaching 
their individuai end-points (convulsion or serious symptom). 

Suppose that n divers are tested simultaneously in a cham- 
ber for a fixed depth and fixed time T, and that the times to 
serious symptom or convulsion are recorded. Since T is a 
fixed time, say 60 minutes, all of the divers will not reach 
their individual end-points (some will complete the entire 
experiment without incident). How does one calculate the 
mean-time-to-failure (mean-time-to-symptom). Let's assume, 
Since the data suggests, that MTTF is exponential. 


Let 
= min([X, ,T] 


where: 


ae 
I 


Fixed time of chamber run or experiment 


ms 
Il 


Time to end-point of ith diver 


The maximum likelihood estimator for MTTS = \ 1s: 


~A[X, 6,42 (1-65) ] 04 A [X6.+T (1-55).] 65 “A [X_ 6 +E (1-9, } 6. 
e A e Mo ue r 


L(,X) 


n ~A[X,6,+T (1-6, )] 6. 
= [If e r 
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where: 


aL 
6. «= 
ale 
0 ede a> AL 
30 
Nn 
Q(X,X) = -A } [X,6.4+T41-6.) + Slog -AX!' + §S log A 
— al 5b ot at ie 
where: 
a ee — log (17 X) 
Nn 
So] 2s 
=a 
Now differentiate: 
oe . ae, S 
oA! eT 
ae | 
3° A 


This leads to: 


hs > 


AS 


Number of failures observed in the interval (0,T) 
The total cumulative exposure time of all divers 


The variance of \ is: 
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var(A) = (E[-S})7* 
r 
= (~ B(s)17* 
r 
= 1 - ~ 
= [5 nfl-e **})™ 
r 
= at 
= ae 
Note that as 
a 2 
ieee fT > ©, Var(XrA) > a which is the same as uncensored. 
Pe \2 1 
ae > 60, eee) x a it >o-efor nf iced. 


wwekt ALT is “small," tcc is small and the precision 
of the maximum likelihood estimator is small. Var(A\) 


is very large, unless n is large. [Ref. 17] 


This complicated mathematical path reduces to a simple ratio 
which is intuitively clear and simple. 


Number of Serious Symptoms Observed in Experiment 


Total Cumulative Exposure Time of All Divers 
[Ref. 16] 


C. CONFIDENCE LIMITS FOR THE EXPONENTIAL PARAMETER 
Suppose that data is available that one wishes to refer to 
an exponential distribution. Let tyetor-- ety, denote the 


observed "times" to symptom or convulsion. Assuming that the 
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times are distributed exponentially one can establish confi- 
dence limits for the rate parameter \X, or for 1/A = E[(T]), 
the mean-time-to-convulsion or symptom. Procedure for two- 


sided confidence limit: 


_ aes 
a Compute t = = ) Ge 
noo. 1 
shel 
b. Decide upon the significance level a. 
eee Conftiaence limits for \’} (failure rate): with confidence 
(l-a~)1003%: 
2 2 
1 *df,a/2 1 *df,1-(a/2) 
—(—=4 J) «k« [eee] 
a elie = ele 


dad. Confidence limits for the Mean-time-to-symptom: 


Z 
~1 — *2n,a/2 — 


2 
X - 
2n,l=(a72) <E(T] = : ee me 


tl Zi 


(Ref. 17] 


D. PARAMETER ESTIMATION 
The probability density of the time to serious symptom 


has been defined.as: 


t 
- f A(t) dt 
f(t) = A(t)e 0 


For the moment let us denote \(t) as a constant failure rate 


y, for a given depth, where y > 0 and substitute y into the 
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above equation. 


Thus we observe that the distribution of failure times is an 
exponential distribution when it can be assumed that the 
failure rate is constant. For this reason, the assumption 
of a constant failure rate is often called the "exponential 
assumption" in reliability models. Now the mean-time-to- 
symptom (MTTS) equals l/y. Of course, in this situation y 1s 
not a constant but is a function of depth, PPO2, PPCO2, work 
rate and perhaps several renee variables. We can however, 
calculate MTTS for those experiments that have been conducted. 
Then with those experiments that have commonality (wet no-work, 
wet working, dry, etc.), project a curve and try to discover 
meaningful relationships from the curve. First, using Donald's 
data since his experiments have the most data points, one ob- 
serves that if the depth is held constant, the time-to-symptom 
is plausibly exponential for any given depth. Moreover, uSing 
the censoring expectation technique we can compute MT'TS Lor 
each experiment. Table 2 tabulates the MTTS for Donald's wet 
resting dives. 

Since y = 1/MTTS, ? can be approximated by 1/MTTS. Figure 
2 is a plot of MTTS as a function of depth and time. Taking 
into account the variability of these averages, a smoothed 


curve of mean-times-to-symptom for non-working divers can be 
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TABLE 2 


fon TIME-TO-SYMPTOM USING THE CENSORING TECHNIQUE FOR RESTING 





DIVERS UNDERWATER AS REPORTED BY DONALD [REF. 10]. « = .05 
, 2 2 
FSW MTTS n 2n *2n(a/2) *2n(1l-a/2) 95% CONFIDENCE LIMIT 
40 199.3 31 62 .675 1.39 134.5 < MTTS < 277.0 
50 wae? | AY 82 .714 ee 53.0 < MTTS < 98.7 
41.5 40 80 714 1.33 29.6 < MTTS < 55.2 
70 32.1 15 30 560 Sy 18.0 < MTTS < 50.4 
55.4 20 40 sila 1.48 33.9 < MTTS < 82.9 
80 20.1 26 52 647 1.43 13.0 < MTTS < 28.7 
90 21.8 55 110 742 ae 16.2 < MTTS < 28.3 
100 9.5 16 32 560 1.57 5.3 < MTTS < 14.9 


Peejeected. Such a projected curve is provided at Figure 3. 
Using this curve one can calculate a conservative estimate of 
what the mean-time-to-symptom is for any depth, and therefore 
calculate a reasonable approximation to y. This composite value 
of y is a function of depth, duration of exposure, PPO2, 


PPCO2, and the other variables, so in some respects y 1S a 


point quantitative value of frit). Now making use of the fact 
migat : 
1e 
Se nase 2 = / £.(t)dt 
0 
t = 
ee ns 7-3 
0 
alin 
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fee) 12s the reliability function of an exponential model. 
Thus from the projected curve we can calculate y for any 


depth based on the projected MTTS suggested by the curve. 


B. MODELS 

If one ponders Figure 3, one can see several approaches 
to model the shape of the curve in the given range. First 
of all it looks exponential or hyperbolic. Now, what we are 
seeking 1S a simple formula in terms of depth and time which 
yields reasonable probabilities for mission success or mission 
failure or the probability that the diver will get a hit. 
One such simple formula that yields a conservative vaiue for 


MTTS in terms of depth is the hyperbolic function: 


MTTS (| peprHy) = C/ (depth - D)) 


where C and D, are constants. Two curves generated by this 


0 
function are shown in Figure 4. This formula allows us to 


express the reliability function of a wet non-working diver 


(a diver in an SDV for example) in terms of depth and time: 


==) 


D 


0 


EXP [ (iE) | Or D > B 


undefined BONE Di< Do 


where D is depth and v is a constant. The relationship between 
such a function and that of a working diver will be explored 


later in this thesis. 
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i Hypenboelic Model l 


Let: 
D = 33FSW 
Vv = ‘laeeZ5 
Gy = 2660 


The following values can be obtained from the formula: 





+. c+? 
(Det DS) T 
FSW MTTS = T Y 
110 8.32 .1200 
100 10.00 .1000 
90 12.39 .0807 
80 16.01 .0625 
70 21.98 .0456 
60 33.37 .0300 
; 50 61.60 .0162 
40 200.00 .0050 
30 ae -- 
Now: 
-(D-D,))” 
RCE 1 fe t] 


= Mission Reliability of a wet non-working 
diver 


ou 





Let: 


t = 20 minutes and depth = 85FSW 
1.325 
_ 65336" 
R(t) = EXP[ = (20) } 
- gamora eo) 
= e 
= .24 


So the probability of mission success under these circumstances 
(say a diver in a wet submersible (SDV) at 85FSW for 20 
minutes) is 24% and accordingly, the probability diver X 


incurs a serious symptom is 76%. 


2. Hyperbolic Model 2 


Eeu.: 
D = 25FSW 
vy = 1.86 
C = 30,800 
FSW MTTS = T d 
100 10.00 LOO” 
90 sec e ey .0765 
80 Sooo OG 
70 25.90 Oe © 
60 Ael..00 .0244 
50 7720 0 Odes 
40 20020 0 . OBO 
30 Spey GS ae .00065 
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Let: 


t = 20 minutes and depth = 85FSW 
iOS 
7 - (85-25) 
R(t) = EXP[ 30,800 (20) | 
_ - .0659 (20) 
= e 
—e 


Mission success is 27% and the probability that diver X gets 
fart 1s 73%. 
3. Exponential Model l 
Another approach to modeling the smoothed projected 


Gmave is to fit an exponential. 


-a (D-Dp) 
T = METS (depth) ~ = ke 
Figure 5 projects such a fit. Let: 

k = 425 

D = 25FSW 

a = .05 

FSW MTTS OL 

100 Das OOH: 
90 16.48 MOO 
80 PA | .0368 
70 44.79 0223 
60 WS ahs NOW 3 5 
50 ie. 1 © 0082 
40 2007 > 20050 
30 SSO SiS .0030 
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Btc) = “ERP ([- =a(D=D)) (t) J 
ke 
Let: 
t = 20 minutes and depth = 85FSW 
R(t) = EXP (-—————Sergemyey (20) | 
425e ° 
aia? oS 20) 
= 
ES 


So the probability of mission success under these conditions 

is 39% and the probability diver X gets a hit is 61%. Given 
the dangers involved with oxygen diving, Hyperbolic Model 1 
provides the most conservative estimates in favor of diver X, 
given Donald's data. Figure 7 shows both the exponential fit 
and the hyperbolic models with the mean-time-to-serious-symptom 


as reported by Donald. 


F. VARIABLE DEPTH DIVE PROFILE 

If we make the exponential assumption about the distribu- 
tion of the failure times (times-to-serious-symptom), some 
very useful results can be derived concerning the mean-time- 
to-symptom (MTTS). Making use of the fact: 


ie 
R(t) = L-F,(t) = L- J £,,(t) dt 
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we obtain 


e 
Re) cel | ye isae =e @ <* 
0 


for the reliability function of the exponential model. 

Suppose that a dive profile consists of n component dives 
at various depths. Consider these dives as connected in a 
series and that the component dives have respective failure 
rates of Vqrlgr +927, If the n component dives are thought 
of as a series, namely the performance of any one component 
dive does not affect the probability of success of any other 
component dive, then under these conditions the probability 
that the dive profile is a success will be the product of the 
component dive reliabilities. 


n 
R IT R. [Ref. 18] 


(PROFILE) 
For example, suppose a dive profile consists of 200 

Minutes at 40FSW with an excursion to 60FSW for 15 minutes. 

What is the probability of success of the dive profile? Using 


model one to determine the value of the y's, 


1 326 1.325 
Uiesay (50288) 
R (PROFILE) EXP [-~——3 639 (200) J EXP | 7630 19) 1 
BGGOR. 8630 
a vy, 
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The probability of mission success with this dive profile is 
2.7% and the probability that diver X gets a hit is 57.3%. 
These are not very good odds for the SDV pilot who must make 
Such an excursion. 
Rt Stewenncce Ormver~ Gn these@xygen Tolerance Curve 

It has been generally accepted that work diminishes 
tolerance to oxygen at increased partial pressures. This was 
demonstrated by Donald [Ref. 10] in the first large series of 
experiments on human beings. His subjects worked vigorously 
underwater by lifting a large bag of weights by means of a 
pulley. He conducted experiments at 50, 40, 35, and 25 FSW. 
These experiments show conclusively that oxygen tolerance is 
markedly diminished by work. This author was unfortunately 
Only able to acquire the data from the experiments at 50 and 
40 FSW. The corresponding-means are 41.5 minutes for 40 FSW 
and 21.9 for 50 FSW. Although this provides only two data 
points, one could hypothesize from the previous analysis and 
project a curve based on the projected curve of the resting 
divers, and fit a model similar to the models proposed earlier 
in this thesis. Figure 7 is an attempt to develop such a 
model. Consider the curve to be asymptotic to 20 FSW. Then 


we have: 
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FSW PROJECTED. CURVE ElteeD CURVE 


20 --- = 
30 eo 0 160.0 
40 Aa 42.0 
50 ie 0 OS 
60 12..0 {ee 8 
rae) ae dae 
80 De 0 Die 
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100 S220 D9 





Figure 8 is a plot of the old and new projected curves with 


the fitted curve based on 


T(depth) = C/ (depth-D,) 


where: 


™ 
II 

‘aaa 
eS 
™= 

on 
un 
Un 


oO 
iI 

N 
© 
2) 
Y) 
= 


Let's compare the results of this model with that for a non- 


working diver. Let: 


D = 20 FSW 

Va = 1293 

Cam 137055 

t = 20 minutes 
depth = 85 FSW 
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(g5-20)1°9? 


cial 13,655 — 


1D, 2) (20) J 


ea 43l(20) 


.0098 


So the probability of mission success under these conditions 
is .98% and the probability the diver X gets a hit is about 
99%. This compares to Res (20) = 24% for the resting diver. 
To this point the analysis and subsequently the pre- 
dictive models have been based on Donald's experiments con- 
ducted in the early 1940's. This is due to the fact that 
Donald conducted the largest series of human experiments to 
date. Since that time the Navy Experimental Diving Unit has 
conducted human experiments but nothing on the scale of those 
conducted by Donald. In the main, experiments since Donald 
have been attempts to probe the limits of the tolerance curve 
of working divers and as seen from Table 3 have not produced 
results similar to Donaid. Most MTTS values are considerably 
larger indicating the tolerance curve is not so restrictive 
as those we have developed thus far. The question arises, 
"What could account for such a disparity?" One possibility 
to be considered is that the breathing apparatus used by Donald 
was not as efficient as those used in later experiments and 
accordingly the PPO2, PPCO2, PPN2 and breathing resistance 
were not equivalent to those of later experiments. Higher 


levels of CO2, N2 and more breathing resistance might account 
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TABLE 3 


MeaN—-TIME=TO=SYMPTOM USING THE CENSORING TECHNIQUE FOR 
DIVERS UNDERWATER PERFORMING WORK 


as 


FSW | EXPERIMENTER MTTS AVERAGE MOST RESTRICTIVE 
20 BUTLER 1016 1447 1016 
1850 
1475 
25 BUTLER 471 TSE0.2 405 
jee 
LANPHIER 405 
PIANTADOSI 1318 
975 
30 BUTLER SS22 Manes no hits 
YARBROUGH 2151 2458.5 aoe 
2007 
LANPHIER 627 
35 BUTLER 470.1 254.95  20).8 
LANPHIER 39.8 
40 BUTLER 593 Da 8 41.5 
83.6 
72.0 
YARBROUGH 294 
447.67 
LANPHIER 88.25 
DONALD 41.5 
45 LANPHIER a7 8 BS 37.5 
50 LANPHIER 29.6 3055.2 21.9 
BUTLER 864 
DONALD Pl ae 
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for this disparity. This hints at the need for human experi- 
ments with elevated levels of CO2, N2, increased breathing 
resistance, and other introduced variables to isolate the 
major factors. The huge variation in MTTS across experimenters 
indicate that some unstudied variable is indeed very important. 
Moreover, Donald's experiments are consistent as a group 
implying the type apparatus and its associated effects, 
temperature, CO2, N2 or metabolic rate as a function of work 
rate and temperature are the prime candidates for study. At 
face value the major variables appear to be PPO2, PPCO2, PPN2, 
and metabolic rate, where metabolic rate is a function of 

work rate, apprehension, water temperature and other varia- 
bles. Such factors as breathing resistance, age of the diver, 
body fat, blood factors, urine factors, exact time to symptom, 
and seriousness of symptom should be recorded. Ideally all 
dives should be carried to a decisive end-point. 

Table 4 is an attempt to find the asymptote for the 
tolerance curve for working divers (50 watts). Many of the 
dives conducted by Butler were variable profile dives, many 
of which involved long dives at 20 FSW prior to an excursion 
to a deeper depth. If one aggregates these dives (using the 
initial part of the variable depth profile) an MTTS = 2641.88 
minutes can be derived. This is well beyond the limits of 
pulmonary toxicity but serves as an ideal anchor for projecting 
a curve and gives credibility to the assumption that the 
iii aie is approximately 20 FSW. It follows that y = .0003785 


and accordingly the probability of a hit at 20 FSW is 
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TABLE 4 
AGGREGATED EXPOSURES AT 20 FSW FOR EXPERIMENTS 
CARRAGED OU Be BUPLEK Al NEDU 


DIVERS TIME OF DIVE TEM EOP A TT TOTAL TIME NOMBERSZOF HLTS 


——_ mm gg al gg mm 


9 x 20° min t 5 = 2285 1 
im x 240 + 304147+45+236 = 3098 * 
16 x 20 = 1920 0 
iLO x 20 = 2280 0 
16 x 240 = 3840 0 
18 x 120 = 2160 0 
ae x 240 Tea eset tO = 9 oie S 

Zio 8 


The aggregated total exposure at 20 FSW equals 21,135 minutes 


during which 8 toxic episodes occurred. 


Beencoreaae | ee 





e7 9003785 (t) 


1- Rog (t) 
abuts) PROBABILITY OF A HIT AT 20 FSW 
60 22% 
120 4.4% 
180 6.5% 
240 8.6% 
300 ORs 
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Hit in this case does not mean convulsion but the symptoms 
as reported by Butler as probable and definite toxic epi- 
meees. Although not precursors of a convulsion these symptoms 
are indicative of the imbalance in the oxygen processes of 
the human body induced by the changes in atmosphere. A 
conservative curve would certainly pass through this point. 
Let us develop one more model using the best guess 
hypothesis. Postulate that the diver's apparatus (SCUBA) is 
working efficiently and that divers control their work rate 
perfectly (a smooth even work rate of about 50 watts). From 
Table 3 our best guess of what the divers could reasonably 


attain might be as follows: 


FSW OBTAINABLE MTTS PIESED GERVE 
20 2641 =e 
30 1000 SBS te Bees, 
40 So0 8270 
0 150 Iho eyae 
60 85 oi es 
rt 50 94.43 
80 5 34.82 
90 25 23.44 
100 le, 16.45 


Figure 9 is a plot of these new curves where DEPTH = slope x MTTS 


to the power of m, 1.ée., 
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D = b(mrrs)™ 


has been fitted to the projected curve. It follows that, 


MTTS = (D/b)2/™ 

where: 

b = 230.16 

l/m = -3.36 

Now 

R(t) = EXP(-——,=(t) | 

(D/b) 
= EXP| (By 7/™ ty] 


Let t = 20 minutes and depth = 85 FSW. 


20) We 8) 


EXP [- (£20258) ~ +39 (20) |] 


Req (20) 


e7 93519 (20) 


~4946 


Thus the probability of mission success is about 49% and the 
probability that diver X gets a hit is about 51%. Let t = 60 


minutes at a depth of 30 FSW. Then 
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230.16)-3 
30 


36 
Ro, (60) EXP [- ( Poou 


=. OOnmew one’) 
e 
= 36 


Here we have the probability of a successful mission for one 
hour at 30 FSW at 93.8%. This curve is for a well trained 
diver whose equipment is working ae His chances of 
a hit is only about 6.2% including symptoms which are not 
abortive. 


Let's examine 120 minutes at 30 FSW. 


230.16)—-3.36 
30 


R3, (120) Jee) || i CLEZ Or} 


7° 90106(120) 


= .88 
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fit. DiSCUSSsaON 


Observation is the basis of our understanding of the world 
around us. But observation only provides information about 
the specific events which we observe: alone, it provides 
little help for dealing with new situations. Useful knowledge 
results from our ability to recognize similarities in dif- 
ferent events, isolate the important factors, and generalize 
from our experience. Roughly speaking, a model is an analogy 
for some real-world phenomenon or process we wish to charac- 
terize. Based on our observations a mathematical model charac- 
terizes the phenomenon in a mathematical structure which 
ideally represents everything that iS important about the 
phenomenon or process and ignores everything else. It is 
difficult to realize this ideal because often we are not sure 
what aspects of the real-world are important. Our goal has 
been to provide a basis for a reasonable and fruitful corres- 
pondence between the observed phenomenon of oxygen toxicity 
in man and the mathematical structure which we have used to 
represent the phenomenon. The objective of this study was 
to integrate the known data points to generate oxygen tolerance 
curves in terms of depth-time limitations. This was done by 
applying a known mathematical model and extrapolating to new 
Geneitiens. TOovaeecomplish this, certain simplifying assump- 


tions were made to make the mathematics tractable and 
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understandable. For the sake of clarity these assumptions 


are listed below. 


A. 


il. 


ASSUMPTIONS 


There iS a distribution of times-to-serious-symptom. 
This distribution looks gamma with a conspicuous 
build-up of hits followed by a long drawn-out tail. 

If one considers the time required to press the diver 
to depth and truncates this part of the empirical 
distribution curve, then the times-to-serious-symptom 
is plausibly exponential. Therefore, the exponential 
distribution provides a reasonable fit to the failure 
Pine TdisGimibwicion fon-a fixed depth. 

Given that the exponential distribution is assumed as 
the distribution of times to serious symptom, then for 
a given depth the parameter for this distribution can 
be estimated by 1/MTTS = - which is a composite value 
of the effects exerted by all the variables important 
GImenemaistribution. By computing MTTS for the experi- 
ments which have been performed, and grouping those 
experiments for which commonality can be assumed, one 
can project a curve from which a value of : can be 
predicted for any depth. 

We have assumed that the phenomenon of oxygen toxicity 
is well behaved allowing us to hypothesize about the 


shape of the curve between data points. This 


Sf 





assumption allowed us to project the curves used in 

this study. 

We have assumed that the integration of the distribution 
of times-to-serious-symptom for a given depth is too 
complicated a process for a model based on the given 
data, and that we could develop a simple formula in 
terms of depth and time which is equivalent to this 
integration, the value of which is a composite value 

of all the variables expressed in the curve. Sucha 
formula provides a reasonable approximation to the value 
that would be obtained from the integration of this 
unknown distribution. 

We have assumed that the major factors in the oxygen 
toxicity phenomenon are depth, time and work. The 

major factors are probably PPO2, time and metabolic 
rate. However, past experiments have not provided 

data precise enough to derive these factors from depth 
and work rate. Moreover, values for PPCO2, PPN2, and 
work rate were unfortunately not always recorded. At 
the very least, future experimenters should attempt 

to measure and record these variables. 

By serious symptom we mean a toxic episode which would 
lead to the abortion of a dive. Many of the symptoms 
used in this study may not have been this serious. 
However, not knowing the degree of impairment associated 


with the reported hit, this author chose to be 
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conservative and use the symptom that ended the experi- 
ment as a serious symptom. Being a Naval Special War- 
fare diver the author is fully aware of the conserva- 
tive nature of this approach. When the dive is long 
and cold, and many dives have already been performed, 
sometimes the fastest way to the bar is not the best 
way for the experimenter or statistician. This builds 
a conservative factor into the models automatically. 

We have assumed that a dive profile can be thought of 
and modeled as a dive to a specified depth followed 
instantaneously by a dive to another specified depth 
and so on, thereby suggesting that a dive profile can 
be modeled as a series of connected, but independent 
dives. This assumption may be hard for some to accept, 
namely that the performance of any one component dive 
does not affect the probability of success of any other 
component dive. The feeling is that there must be some 
sort of cumulative oxygen dose that one builds up. 
However, if one assumes a series of dives, as we Nave, 
there is indeed a cumulative effect in that the proba- 
bility of a successful dive profile is the product of 
the individual component dive probabilities. 

Individual variation can be described by the normal 
distribution. Statistical theoreticians have noticed 
that the variance associated with complex behavioral and 


biological phenomena is often approximately normally 


ae, 





distributed. Moreover, any time a phenomenon is depen- 
dent upon a number of underlying factors, the phenomenon 
itself will be normally distributed, regardless of the 
shape of the underlying distributions. This appears 

to be the case in oxygen toxicity. Donald tested the 
same man on 20 separate occasions and found a huge 
variation in the time-to-symptom. What the underlying 
factors are remain to be discovered and are probably 
the key to developing an accurate tolerance curve. 
However, if one assumes that there is a time about 
which there are equally distributed, the same number 

of hits for a given depth, then normality is plausible, 
and indeed more than likely. What we see in the 
aggregate of individuals does not look normal but gamma. 
We can still assume normality for individual variation 
and exponentiality for the aggregation of all divers 
because some individuals appear to be more toierant 

than others to higher than normal partial pressures of 


oxygen. 


Be RESULTS 

We were able, given the above assumptions, to mathematically 
model the reported incidence of oxygen toxicity. The proba- 
bility of mission success, which here is defined as the proba- 
bility the diver can accomplish the dive profile without a 


toxic episode of such a degree as to abort the dive, was 
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modeled as 


c 
= gente) Kets 
R(t) = e 


where f(t) is the distribution of times to symptom. The 
models developed herein were not cross validated: testing 
and further development of this approach must await further 
studies in oxygen toxicity. 

Many investigators will feel uncomfortable with this 
statistical black-box approach to oxygen toxicity, and it is 
understandable that one desires a precise physical concept as 
an anchor. These preliminary results merely show that the 
use of a mathematical model may provide a means to deal with 
the variance associated with oxygen toxicity. Such models 
transcend the tested profiles and should improve the preci- 
sion of estimating toxic episodes and allow the associated 
risk to be evaluated. 

It is stressed that the applicability of probability theory 
to oxygen toxicity has not herein been proven, but a serious 
study based on animals treated in such a fashion could demon- 
strate the usefulness of such models and could provide some 
useful parallels. Moreover, if future researchers test for 
and record all the major variables, statistical analysis may 
provide more accurate models incorporating water temperature, 
work rate, PPN2, PPO2, PPCO2, and human factors such as age, 


percent of body fat, experience, etc. Until such time as the 
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physical concept is discovered or an accurate mathematical 
model is developed, we will continue to operate on trial and 
error tables which are rules analogous to the old 30-feet 
for 30-minute rule. Moreover, a few chance hits, due to the 
variability involved, may preclude the use of profiles which 
are useful, to the prescribed degree of safety, and which are 


needed by the fleet. 
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EV. “CONCEUSTIONS 


SCemeral Nervous System (CMS) oxygen toxicity in man ais 
characterized by wide individual variation and diverse symptoms. 
And although the deleterious effect of oxygen under pressure 
has been recognized for almost a century, the mechanism(s) by 
which high pressure oxygen produces convulsions or other toxic 
manifestations remains unresolved. Several theories have been 
Suggested to explain the cause of CNS oxygen toxicity, but 
no Single hypothesis has received universal acceptance. This 
study attempts to determine the shape of the tolerance curve 
by statistical analysis of the existing data base. 

Statistical theoreticians have noticed that the variance 
associated with complex behavioral and biological phenomena 
is often approximately normally distributed. Moreover, any 
time phenomena are dependent apen = number of underlying fac- 
tors (which appears to be the case with oxygen toxicity). 
The phenomenon itself will be normally distributed, regardless 
of the shape of the underlying distributions. With this 
normality assumption applied to individual variability, and 
the time-to-serious-symptom assumed to be exponentially dis- 
tributed, this study develops oxygen tolerance curves in 
terms of depth-time limitations. 

By applying a known distribution to the time-to-serious- 


symptom, mathematically predictive models have been developed 
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which allow a greater degree of predictability in mission 
profiles while concomitantly allowing the associated risk to 
be evaluated in terms of probability. This is one of the 


major objectives of mathematically modeling any phenomenon. 
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